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ABSTRACT 

Hydrodynamic simulations of growth of cosmic structure suggest that 30-50% of the total baryons 
at z = may be in a warm-hot intergalactic medium (WHIM) with temperatures ~ 10^ — lO^K. The 
O VI AA1032, 1038 absorption line doublet in the FUV portion of QSO spectra provides an important 
probe of this gas. Utilizing recent hydrodynamic simulations, it is found that there should be ~ 5 O VI 
absorption lines per unit redshift with equivalent widths > 35 mA, decreasing rapidly to ^ 0.5 per unit 
redshift at > 350 mA. About 10% of the total baryonic matter or 20 — 30% of the WHIM is expected to 
be in the O VI absorption line systems with equivalent width > 20 mA; the remaining WHIM gas may 
be too hot or have too low metallicity to be detected in O VI. We find that the simulation results agree 
well with observations with regard to the line abundance and total mass contained in these systems. 
Some of the O VI systems are collisionally ionized and some are photoionized, but most of the mass is in 
the collisionally ionized systems. We show that the gas that produces the O VI absorption lines does not 
reside in virialized regions such as galaxies, groups, or clusters of galaxies, but rather has an overdensity 
of 10 — 40 times the average density. These regions form a somewhat connected network of filaments. 
The typical metallicity of these regions is 0.1 — 0.3 Z©. 

Subject headings: Cosmology: large-scale structure of Universe - cosmology: theory - intergalactic 
medium - quasars: absorption lines - numerical method 



1. INTRODUCTION 

Three independent observations at z > 2, namely the 
latest cosmic microwave background experiments (Netter- 
field et al. 2001) at z >> 2, the Lyman alpha forest decre- 
ment measurements (Ranch et al. 1997; Weinberg et al. 
1997) at z ~ 2 — 3, and the deuterium-to-hydrogen ra- 
tio measurements (interpreted within the context of the 
standard nucleosynthesis theory. Buries & Tytler 1998), 
aU consistently give Vlb{z > 2) = (0.019 ± 0.002)^,-2 ^ 
0.039 ± 0.004 (where h = Hq/IOO km sec"! Mpc'^ = 0.7 
is used for the last equality). On the other hand, in our 
local universe all the well-measured baryonic components 
do not appear to add up to the indicated baryonic den- 
sity seen at high redshift by a factor of about three (e.g., 
Fukugita, Hogan, & Peebles 1998) with f]f,(z = 0) = 
+ Q.HI + + ^XTay,ci ~ 0.0068 < 0.011 (2(7 limit), 
where Q.hi, and ^xray,ci are the baryonic densi- 
ties contained in stars, neutral atomic hydrogen detected 
in 21cm or damped Lyman alpha systems, molecular hy- 
drogen and hot X-ray emitting gas in rich cluster cen- 
ters. While a substantial amount of baryons may be in 
the low-redshift Lya forest with f^Lya = 0.008 (Penton, 
ShuU, & Stocke 2000), it is still far short of accounting 
for the baryons seen at high redshift. Thus, most of the 
baryons are "missing" today. 

Hydrodynamic simulations of structure growth suggest 
that a large fraction of the missing baryons at the present 
epoch may be in a gaseous phase with temperatures of ~ 
10^ — 10'' K at moderate overdensities, typically ~ 10 — 40 
(Cen et al. 1995; Cen & Ostriker 1999a; Dave et al. 2001). 
It is important to detect this gas to enable a better under- 
standing of the baryon distribution in the local universe 



and because it may have an important influence on the 
evolution of galaxies (Blanton et al. 2000). The O VI 
AA1032, 1038 absorption line doublet in the spectra of low- 
redshift QSOs provides a valuable probe of the quantity 
and properties of this gas, and recent observations with 
the Space Telescope Imaging Spectrograph on the Hubble 
Space Telescope have shown that the number of interven- 
ing O VI absorbers per unit redshift {dN/dz) is quite high 
(Tripp, Savage, & Jenkins 2000; Tripp & Savage 2000), 
suggesting that these absorbers do indeed trace a signifi- 
cant baryon repository. In this Letter we use cosmological 
hydrodynamic simulations to make detailed predictions of 
the O VI absorption line properties, and we compare these 
predictions to recent observations. This complements ear- 
lier work by Hellsten, Gnedin and Miralda-Escude (1998), 
which focused on higher ionization oxygen lines observable 
with X-ray telescopes. 

2. SIMULATION 

We use two recent cosmological hydrodynamic simula- 
tions of the canonical, cosmological constant dominated 
cold dark matter model (Ostriker & Steinhardt 1995) 
with the following parameters: Qm ~ 0.3, Q\ ~ 0.7, 
rjfc/i^ — 0.017, h — 0.67, as = 0.9, and the spectral index of 
the primordial mass power spectrum n = 1.0. The higher 
resolution simulation box has a size of 25h~^ comoving 
Mpc on a uniform mesh with 768'^ cells and 384^^ dark 
matter particles. The comoving cell size is 32.6h~^ kpc. 
The mass of each dark matter particle is 2.03 x 10^ M0, 
and the mean baryonic mass per cell is 3.35 x 10^h~^MQ. 
An earlier simulation with a coarser resolution but a larger 
box of size L = lOO/i^^Mpc (Cen & Ostriker 1999a,b) is 
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also used to calibrate the boxsize effect. In velocity units 
the cell sizes are A7; = (3.3, 19.5) km/s in the two simula- 
tions, respectively, providing adequate spectral resolution 
in comparison to current observations.^ 

A word on simulation resolution is relevant here. As sim- 
ulations with much higher spatial resolutions have shown 
(Dave et.al. 2001), the bulk of the warm-hot intergalac- 
tic medium has an overdensity of 10 — 40. This has an 
important implications: the WHIM gas is mostly in un- 
coUapsed regions, for which the Jeans mass is well resolved 
by our high resolution simulation. In fact, even for a 
lower temperature, photoionized gas at T = lO"' K, the 
Jeans mass is ^ 10® and Jeans length is 400 kpc/h 
comoving, which are to be compared to the 3.3 x 10^ Mq 
and 2 x lO'^M© nominal baryonic and dark matter mass 
resolution, and 33kpc/h comoving nominal spatial resolu- 
tion of our 25Mpc/h box; the corresponding numbers for 
our lOOMpc/h box are 7.1 x W'^Mq and 5.1 x lO^Mg, 
and 195kpc/h. The true spatial and mass resolutions of 
the simulation are probably a factor of 2 and 8 times 
worse (Ccn & Ostriker 1999b). Clearly, our higher res- 
olution simulation is adequate for resolving unvirialized 
intergalactic gas that was initially photoionized and then 
shock heated to the indicated temperature. Our lower res- 
olution simulation is marginal at resolving smaller WHIM 
structures, which will be evident in the results shown be- 
low. 

We follow star formation and feedback using a well de- 
fined prescription similar to that of other investigators. A 
stellar particle of mass = c*mgasAt/i* is created, if 
the gas in a cell at any time meets the following three 
conditions simultaneously: (i) flow contracting, (ii) cool- 
ing time less than dynamic time, and (iii) Jeans unstable, 
where At is the timcstcp, t^, — max(tdyn7 lO'^yrs), tdyn is 
the dynamical time of the cell, mgas is the baryonic gas 
mass in the cell and c* = 0.25 is star formation efficiency. 
Each stellar particle has a number of other attributes at 
birth, including formation time and initial gas metallicity. 
The typical mass of a stellar particle is one million solar 
masses. Stellar particles are subsequently treated dynam- 
ically as collisionless particles, except that feedback from 
star formation is allowed in three forms: UV ionizing field, 
supernova kinetic energy, and metal enriched gas, all be- 
ing proportional to the local star formation rate. Metals 
are ejected into the local gas cells where stellar particles 
are located using a yield Y = 0.02 (Arnett 1996) and are 
followed as a separate species adopting the standard solar 
composition. The simulation also includes cooling/heating 
processes due to the self-consistently produced metals us- 
ing a code based on the Raymond-Smith code assuming 
ionization equilibrium (Cen et al. 1995). We find that the 
computed metallicity distributions over a wide range of en- 
vironments, including clusters of galaxies, damped Lyman 
systems, Lyman alpha forest and stars, are in broad agree- 
ment with observations (Cen & Ostriker 1999c; Nagamine 
et al. 2001; Cen et al. 2001), lending us confidence that 
the computed metal distribution in the intermediate re- 
gions under consideration (between Lyman alpha forest 
and clusters of galaxies) may be a good approximation to 
the real universe. 

Post-simulation the photoionization code CLOUDY 



(Ferland et al. 1998) is used to compute the abundance of 
O VI, adopting the shape of the UV background calculated 
by Haardt & Madau (1996) normalized by the intensity 
at 1 Ryd determined by ShuU et al. (1999) and assum- 
ing ionization equilibrium. We generate synthetic absorp- 
tion spectra using a code similar to that used to generate 
the synthetic Lyman alpha forest spectra in earlier papers 
(Cen et al. 1994; Miralda-Escude et al. 1996). To eluci- 
date the physical nature of the absorbers, we distinguish 
between photoionized and collisionally ionized gas; O VI 
absorption lines with optical depth-weighted temperature 
less than 10^ K are identified as photoionized systems and 
others as collisionally ionized lines. However, in some cases 
both ionization mechanisms are important. 




^The STIS observations reported by Tripp et al. 2000 and Tripp & Savage 2000 have spectral resolution of 7km/s (FWHM). 
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Fig. 1. — Three-dimensional distribution of the O VI density for 
photoionized gas (top panel) and coUisionally ionized gas (bottom 
panel). The box size is 25/i~^Mpc. The color scheme is as follows: 
light green, green, bright green, green-yellow, yellow, ycUow-rcd and 
red-black regions have O VI densities of lO"''', IQ-'^-^, 10"'', IQ-^-^, 
10-'^-^ and > IQ-^-^ in units of Pcr«- 



3. RESULTS 

Figure 1 shows the 3-d distribution of the O VI den- 
sity: Figure la shows the photoionized O VI and Figure 
lb shows the coUisionally ionized O VI. It is seen that 
the low density (bright green) O VI gas is somewhat con- 
nected to form a filamentary structure, while higher den- 
sity gas surrounds roundish regions. Galaxies usually sit 
in the middle of the roimdish red regions (Cen & Ostriker 
1999c). The O VI lines due to photoionized gas tend to 
stem from lower density regions. As we will show later, 
the bright-green and green-yellow regions produce most 
of the observed OVl absorption lines; these regions have 
overdensities of 10 — 40 times the average density. 
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Fig. 2. — Cumulative number of OVl absorption lines per unit 
redshift as a function of the O VI absorption line equivalent width for 
all lines (solid curve), lines due to photoionized gas (dotted curve) 
and lines due to coUisionally ionized gas (dashed curve) with Icr 
statistical crrorbars. The left panel is from the L = 25/i~^Mpc box 
and the right panel is from the L = 100/i^^Mpc box. The equiva- 
lent widths arc computed by integrating the flux spectrum between 
downcrossing and upcrossing points at 80% of the continuum. Both 
dotted and dashed curves are shifted downward by 5% for clarity. 
The symbols are observational data from Tripp, Savage, &c Jenk- 
ins (2000; solid dot) and Savage et al. (2001; open dot) with Icr 
errorbars. 

Figure 2 shows cumulative number of lines per unit red- 
shift as a function of the O VI absorption line equivalent 

width for all lines (solid curve), photoionized lines (dotted 
curve) and coUisionally ionized lines (dashed curve) , com- 
pared to observations from Tripp et al. (2000) and Savage 
et al. (2001). We note that at EW < 35mA the pho- 
toionized lines outnumber the coUisionally ionized lines. 
But at the high end of EW the coUisionally ionized lines 
dominate in total numbers of lines, due to the fact that 
high EW lines originate only in high temperature and 
higher density regions. Visual examinations of synthetic 



O VI absorption spectra indicate that O VI absorption 
lines due to photo-ionized gas are typically narrower (6 
parameter) than O VI absorption lines due to collision- 
ally ionized gas. The agreement between simulations and 
observations is reasonably good, given the small observa- 
tional samples and computational issues which may have 
caused both boxes to underestimate the number of lines 
at various equivalent widths. The larger simulation box 
(right panel) definitely underestimates the number of lines 
at EW < 70mA due to inadequate resolution, while even 
the smaller box may be underestimating the number of 
lines at EW < 20mA and still higher spatial resolution 
(in progress) would be required to adequately model these 
lines. But, at larger width, cosmic variance is substan- 
tial and the smaller box underestimates the abundance 
of lines with EW > 350mA. We argue that a box of 
size L > lOOMpc/h may be required, if one is interested 
in studying lines with EW > 350mA. We expect that 
improvement of numerical simulations will likely increase 
the computed number of lines at both ends of the equiv- 
alent width distribution. Recently, Fang & Bryan (2001) 
have independently obtained very similar results from a 
cosmological simulation which employs different numeri- 
cal methods. 
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Fig. 3. — Cumulative distributions of the fraction of baryonic gas 
contained in OVl absorption line, again, for all lines (solid curve), 
lines due to photoionized gas (dotted curve) and lines due to col- 
lisionally ionized gas (dashed curve). The left panel is from the 
L = 25/j^"'Mpc box and the right panel is from the L = 100h~^yipc 
box. Both dotted and dashed curves are, again, shifted downward 
by 5% for clarity. The symbols are observational data from Tripp, 
Savage, & Jenkins (2000; solid dot) and Savage ct al. (2001; open 
dot) with Icr crrorbars, assuming {Z/ Zr.))f{OV I) = 0.02 adopted 
from simulation (see the bottom panel of Figure 4) 

Figure 3 shows the cumulative fraction of baryonic gas 
density contained in the O VI lines as a function of 
equivalent width for all lines (solid curve), photoionized 
lines (dotted curve) and coUisionally ionized lines (dashed 
curve), compared to observational constraints. Appar- 
ently, most of the mass probed by O VI lines is in lines 
that are coUisionally ionized, although the photoionized 
lines are dominant in number. This is because collision- 
ally ionized lines originate in higher density regions (see 
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Figure 4 below). We see that the baryonic gas contained 
in O VI hnes is about 10% of total baryonic matter or 
roughly 20 — 30% of the WHIM. The agreement with ob- 
servations is good, together with the agreement found in 
Figure 2 strongly suggesting the existence of the WHIM 
in the real universe and that the inventory of mass derived 
from the O VI absorption lines accounts for a significant 
fraction of WHIM. The remaining 70-80% of the WHIM 
is likely gas which it too hot or does not contain enough 
oxygen to be detected via O VI absorption. 

Finally, Figure 4 reveals the physical properties of the 
underlying gas that produces the O VI linos. Wo see that 
at EW = (35, 100)mA the typical mctallicity of the gas 
that produces O VI absorption lines is 0.1 — 0.3 Zq for 
both photoionized lines and coUisionally ionized lines. The 
typical overdensity is ~ (10 — 40) for both photoionized 
coUisionally ionized lines in the range EW = 20 - 200mA. 
Photoionized lines tend to reside at the lower end of the 
overdensity range at the lower EW end, but two types con- 
verge to a typical overdensity of roughly 40 at the high EW 
end under consideration. Evidently, the gas that causes 
the O VI absorption lines is relatively distant from virial- 
ized regions such as galaxies, groups, or clusters of galaxies 
in truly intergalactic, uncollapsed regions. 
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Fk;. 4. — The physical properties of the underlying gas that 
produces the O VI lines. The top panel shows mctallicity of O VI 
lines as a function of equivalent width for lines due to photoionized 
gas (dotted curve) and lines coUisionally ionized gas (dashed curve), 
with Icr dispersion. The symbols arc observed metallicitics derived 
by Tripp et al. (2001) and Savage et al. (2001) assuming the gas 
is photoionized (solid dots) or coUisionally ionized (open dot) with 
2(T errorbars. The middle panel shows the overdensity of regions 
that produce the O VI absorption lines as a function of equivalent 
width for lines due to photoionized gas (dotted curve) and due to 
coUisionally ionized gas (dashed curve) with Icr dispersion. The bot- 
tom panel shows the product of metallicity and O VI fraction as a 
function of equivalent width for lines due to photoionized gas (dot- 
ted curve) and lines coUisionally ionized gas (dashed curve) with Icr 
dispersion. This quantity is useful for calculating Qf,{OVI) from 
observations, see eqn. 1 of Tripp et al. (2000). 



4. DISCUSSION AND CONCLUSIONS 

We have made detailed analyses of our latest hydrody- 
namic simulations to examine the expected properties of 
O VI absorption lines due to intergalactic gas. In the sim- 
ulations the vast majority of O VI absorption lines with 
equivalent width of > 35mA are not due to gas that resides 
in virialized regions such as galaxies, groups, or clusters of 
galaxies, but rather due to intergalactic gas of an overden- 
sity of 10 — 40. These regions form a somewhat connected 
network of filaments. The typical metallicity of these re- 
gions is 0.1 — 0.3 Zq. a relevant test for this contribution 
is to check in both the real and simulated sky: what is 
the distance between O VI absorption systems and bright 
galaxies? 

We find that the simulation results agree well with ob- 
servations with regard to the line abundance and total 
mass contained in these systems. At O VI absorption line 
cqtiivalent width of > 35mA one expects to see ~ 5 lines 
per unit redshift. The number of lines decreases rapidly 
towards larger equivalent width, dropping to about 0.5 per 
unit redshift at > 350mA, above which our current simula- 
tions may be significantly underestimating the line abun- 
dance. About 10% of total baryonic matter or 20 — 30% 
of WHIM is expected to be in the O VI absorption line 
systems with equivalent width > 20mA. Clearly, a larger 
observational sample of O VI absorption lines would be 
highly desirable to increase the statistical accuracy as well 
as other independent observations including absorption 
lines due to other species and direct emission measure- 
ments. 

The evolution of WHIM is likely a function of cosmo- 
logical models, and to some extent, parallels the evolution 
of clusters of galaxies thus provides a powerful, indepen- 
dent test of cosmological models. The model appears to 
be qualitatively consistent with the latest observations at 
higher redshifts (Reimers et al. 2001). We will explore this 
issue of evolution of WHIM quantitatively in a separate 
paper. 
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